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ABSTRACT 

We present the first extensive study of the coronal line variability in an active 
galaxy. Our data set for the nearby source NGC 4151 consists of six epochs of quasi- 
simultaneous optical and near-infrared spectroscopy spanning a period of about eight 
years and five epochs of X-ray spectroscopy overlapping in time with it. None of 
the coronal lines showed the variability behaviour observed for the broad emission 
lines and hot dust emission. In general, the coronal lines varied only weakly, if at all. 
Using the optical [Fe vii] and X-ray O vn emission lines we estimate that the coronal 
line gas has a relatively low density of n e ~ 10 3 cm -3 and a relatively high ionisation 
parameter of log U ~ 1. The resultant distance of the coronal line gas from the ionising 
source is about two light years, which puts this region well beyond the hot inner face 
of the obscuring dusty torus. The high ionisation parameter implies that the coronal 
line region is an independent entity rather than part of a continuous gas distribution 
connecting the broad and narrow emission line regions. We present tentative evidence 
for the X-ray heated wind scenario of Pier & Voit. We find that the increased ionising 
radiation that heats the dusty torus also increases the cooling efficiency of the coronal 
line gas, most likely due to a stronger adiabatic expansion. 

Key words: galaxies: Seyfert infrared: galaxies - X-rays: galaxies - quasars: emis¬ 
sion lines - quasars: individual: NGC 4151 


1 INTRODUCTION 

In addition to the broad and narrow emission lines, the spec¬ 
tra of active galactic nuclei (AGN) display high-ionisation 
emission lines, the so-called coronal lines, which require en¬ 
ergies > 100 eV to be excited. The coronal line region 
is believed to he at distances from the central ionising 
source intermediate between those of the broad (BELR) 
and narrow emission line region (NELR) and to possi¬ 
bly coincide with the hot inner face of the cir cumnuclear, 
obscu ring dusty torus (as first suggested by IPier fe Voitl 
1199511 . Support for this assumption comes from the fact 
that the coronal lines have higher critical densities for col- 
lisional deexcitation than the low-ionisation narrow emis¬ 
sion lines (n e ~ 10 7 — 10 10 cm -3 ), their profiles tend 
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to have full widths at half maxima (FWHM) interme¬ 
diate between those of the broad and na rrow emission 
lines (FWHM~ 500 — 1500 km s -1 ; e.g.. IPenston et alj 
198jj_ Aggenzcllcr & Ocstrciche3_ 1988 : lErkens et al.l Il997l : 


Rodrig uez-Ard ila et al. 2002, 201lj) and the emission from 


this region is often extended but much less so than that 
from the lo w-ionisation NELR (on scales of ~ 80 — 
150 pc; e.g.. iPrie to et al .ll2005l: iM iiller Sa nchez et ai]|2006l : 


iMuller-Sanchez et al. 2011 : iMazzalav et al. 2013h . Coronal 
lines are observed with smilar freq uency in both types of 
AGN llOsterbrockll 19771 : |Koskilll978l l. but type 1 AGN have 
stronger coronal line emission relati ve to their low-ionisation 
narro w lines than type 2 AGN (jMuravama fe Taniguchil 
119981) . Therefore, it is likely that the coronal line region has 
two components, one compact and one spatially extended, 
with only the latter remaining unobscured by the dusty torus 
in type 2 AGN. 


The high ionisation potentials of the coronal lines can 
be produced either in a hot, collisionally ionised plasma, as is 
the case for the solar corona from which these lines have their 
name, or in a gas photoionised by the hard continuum of the 
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AGN. In the first case, the electron temperatures would be 
of the order of T e ss 10 6 K and in the second case much lower 
(T e ~ 10 4 — 10 5 K). Currently, photoionisation is favoured, 
since for most AGN the observed flux ratios between differ¬ 
ent coronal lines c an be reproduced within a factor of ~ 2— 3 
by these models (lOliva et al.l ll9941 : iFerguson et al.l Il997t) , 
whereas the temperature of the hot plas ma wou ld need t o 
be fine-tuned within a very narrow range dOliva et al -1993)- 
In any case, the coronal line region is most likely dust free, 
since strong emission from refractory elements such as iron, 
silicium and calcium are observed, which would be severly 
reduced in a dusty environment. 

The coronal lines are often bl ueshifted relative to 
the low-ionisation narrow lines (e.g.. IPens ton et all Il984l : 
lErkens et al.lll997l : iRodriguez-Ardila et alJl2002h . which in¬ 
dicates that the coron al line gas is in outflow. However, as 
iMullanev et al.1 (i2009l l concluded for the source Ark 564, 
the potential coronal line emitting clouds must have under¬ 
gone acceleration to the observed velocities prior to these 
lines being produced, i.e. they have reached their terminal 
velocity against the opposing drag and gravitational forces. 
Then, given also the similar estimated physical conditions 
and location, it has been proposed that the partly ionised 
gas that produces the O VII and O VIII absorption lines and 
edges seen in the soft X-ray spectra of many AGN, i.e., the 
so-called ’warm absorber’, produces also the coronal lines in 
its (colder) outer r egions IlNetzerlIl993l : lErkens et al.l 1 19971 : 
IPorauet et al.l[l999l l . In this wind model for the coronal line 
region, a considerable contribution from shock ionisation is 
expected. 

The most stringent constraints on the properties of the 
coronal line emitting region could come from variabililty 
studies, in particular if the variability of several coronal 
lines can be compared with each other and with that of 
other AGN components such as the BELR and the X-ray 
continuum. However, mainly due to the weakness of these 
emission lines and also lack of data , very few studi es of this 
kind have been attempted so far. IVeilleuxl dl988l ) did the 
only systematic study of the coronal line variability. In his 
sample of ~ 20 AGN he found firm evidence that both the 
[Fe Vll] A6087 and [Fe x] A6375 emission lines varied (dur¬ 
ing a period of a few years) for only one source (NGC 5548) 
and tentative evidence for another seven sources (includ¬ 
ing NGC 4151). Then, within a general optical variability 
campaign on the source Mrk 110 lasting for half a year, 
iKollatschnv et akl (l200ll ) reported strong [Fe x] variations. 
More recently, follow-up optical spectroscopy of a handful 
of objects with unusually prominent coronal lines selected 
from the Sloan Digital Sky Survey (SDSS) showed that in 
half of them the coronal lines strongly faded (by factors of 
~ 2 — 10), making these sources candidates for stellar tid al 
disruption events llKomossa et al.ll2009l : lYang et al.ll2013fl . 

In this paper, we present the first extensive study of 
the coronal line variability in an AGN. Our data set for 
the nearby, well-known source NGC 4151 (z = 0.0033) 
is unprecedented in that it includes six epochs of quasi- 
simultaneous optical and near-IR spectroscopy spanning a 
period of ~ 8 years and five epochs of X-ray spectroscopy 
overlapping in time with it. Furthermore, the observations 
in each wavelength were performed with the same telescope 
and set-up. This paper is organised as follows. In Section 
2, we present the data and measurements. In Section 3, we 


discuss the observed variability behaviour of the near-IR, 
optical and X-ray coronal lines, for which we seek an inter¬ 
pretation in the context of the location of and excitation 
mechanism for the coronal line emission region in Section 
4. Finally, in Section 5, we summarise our main results and 
present our conclusions. Throughout this paper we have as¬ 
sumed cosmological parameters Hq = 70 km s _1 Mpc” 1 , 
17 m = 0.3, and 17 a = 0.7. 


2 THE DATA AND MEASUREMENTS 


2.1 Near-IR and optical spectroscopy 

We have six epochs of quasi-simultaneous (within 3-14 days) 
near-IR and optical spectroscopy for NGC 4151 (see Tables 
m an d ® ■ The near-IR spectroscop y was obtained with the 
SpeX spectrograph dRavner et al.ll2003l ) at the NASA In¬ 
frared Telescope Facility (IRTF), a 3 m telescope on Mauna 
Kea, Hawai’i, in the short cross-dispersed mode (SXD, 
0.8 —2.4 gm). All data except those from 2010 were obtained 
through a slit of 0.8 x 15" giving an average spectral resolu¬ 
tion of full width at half maximum (FWHM) ~ 400 km s -1 , 
A narrower slit of 0.3 x 15" was used for the 2010 epoch. The 
four epochs spanning t he years 2004 — 20 07 arc ou r own data 
and w ere presented in lLandt et al.l ll2008ll and lLandt et al.l 
ll201lll . The near-IR spectra from 2002 and 2010 were dis¬ 
cussed by iRiffel et~akl (l2006l l and ISchniille et al.l (l2013l l , re¬ 
spectively. Th e optical spectra were obtained with the FAST 
spectrograph dFabricant et al.ll 1998h at the Tillinghast 1.5 m 
telescope on Mt. Hopkins, Arizona, using the 300 1/mm grat¬ 
ing and a 3" long-slit. This set-up resulted in a wavelength 
coverage of ~ 3720 — 7515 A and an average spectral reso¬ 
lution of FWHM ~ 330 km s _1 . The slit was rotated to the 
parallactic angle only for the January 2006 epoch, however, 
except for the 2004 data, all spectra were observed at a very 
low airmass (sec z ~ 1.05). The May 2004 spectrum was ob¬ 
served at an airmass of sec z ~ 1.3 and so the flux loss due to 
atmospheric differential refraction is expected to be ~ 20% 
at the observed wavelength of [Fe Vll] A3759 relative to that 
at wavelengths 5000 A (IFilippenko|[l 98^1. The January 
2006 data were discussed in Landt et all ^ 2008l l. whereas all 
other optical spectra were retrieved from the FAST archive. 

We have measured the fluxes of the strongest near-IR 
and optical coronal lines by integrating the observed pro¬ 
files over the local continuum, i.e. we have not assumed a 
specific line shape. In the near-IR, we have measured two 
sulfur lines and two silicon lines, namely, [S VIIl] A9911, 
[S ix] 1.252/nn, [Si Vi] 1.965/rm and [Si x] 1.430/rm (see 
Table [T]). Two of these emission lines are blended, namely, 
[S ix] with [Fe II] 1.257/rm, and [Si Vi] with H 2 1.957/iin. 
In these cases, we have assumed a Gaussian profile in the 
deblending procedure. In the optical, we have measured four 
iron emission lines, namely, [Fe Vll] A3759, [Fe Vll] A5159, 
[Fe vn] A5721 and [Fe vn] A6087 (see Tablc0. We have con¬ 
sidered also [Fe x] A6375, but this line is extremely weak in 
NGC 4151 and makes up o nly ~ 15% of the total blend with 
[O 1 ] A6364 rtPelat et alJll987u . We do not observe significant 
variability in the flux of the entire blend, which constrains 
the [Fe x] variability to below a factor of ~ 2. 

Since both the near-IR and optical spectra were ob¬ 
tained in non-photometric sky conditions, we study in the 
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Table 1. Near-IR emission line fluxes and ratios 


Observation 

[S III] A9531 

[S III] A9069 

[S III]/ 

[S VIII] A9911 

[S III]/ 

[S IX] 1.252/rm 

[S III]/ 

Date 

(erg/s/cm 2 ) 

(erg/s/cm 2 ) 

[S III] 

(erg/s/cm 2 ) 

[s viii] 

(erg/s/cm 2 ) 

[S IX] 

2002 Apr 23 

(7.61±0.09)e—13 

(2.66±0.02)e—13 

2.86T0.04 

(2.49±0.21)e—14 

30.6T2.6 

(2.44±0.12)e-14 

31.2T1.6 

2004 May 23 

(1.06±0.02)e-12 

(3.89±0.08)e—13 

2.73T0.08 

(5.35±0.60)e—14 

19.8T2.3 

(2.64±0.28)e—14 

40.2T4.3 

2006 Jan 8 

(1.17±0.01)e—12 

(4.40±0.06)e—13 

2.66T0.04 

(3.81±0.26)e—14 

30.7T2.1 

(3.41±0.14)e—14 

34.3T1.4 

2006 Jun 12 

(9.97±0.10)e—13 

(3.85±0.07)e—13 

2.59T0.05 

(3.38±0.19)e—14 

29.5T1.7 

(3.54±0.12)e—14 

28.2T1.0 

2007 Jan 24 

(8.17±0.14)e—13 

(2.79±0.06)e—13 

2.93T0.08 

(3.09±0.38)e—14 

26.4T3.3 

(3.01±0.29)e—14 

27.1T2.6 

2010 Feb 27 

(8.89±0.15)e—13 

(2.98±0.06)e—13 

2.98T0.08 

(4.00±0.36)e—14 

22.2T2.0 

(3.31±0.19)e—14 

26.9T1.6 

Observation 

[S III] A9531 

[Si VI] 1.965/rm 

[S III]/ 

[Si X] 1.430/rm 

[S III]/ 

Pa/3 

[S III]/ 

Date 

(erg/s/cm 2 ) 

(erg/s/cm 2 ) 

[Si VI] 

(erg/s/cm 2 ) 

[Six] 

(erg/s/cm 2 ) 

Pa/3 

2002 Apr 23 

(7.61±0.09)e—13 

(6.75±0.23)e—14 

11.3T0.4 

(4.75±0.22)e-14 

16.0T0.8 

(9.22±0.09)e-13 

0.825T0.013 

2004 May 23 

(1.06±0.02)e—12 

(9.91±0.51)e—14 

10.7T0.6 

(6.23±0.58)e—14 

17.0T1.6 

(7.59±0.24)e-13 

1.397T0.053 

2006 Jan 8 

(1.17±0.01)e-12 

(8.56±0.33)e—14 

13.7T0.5 

(4.51±0.19)e—14 

25.9T1.1 

(l.ll±0.01)e—12 

1.054T0.014 

2006 Jun 12 

(9.97±0.10)e—13 

(7.59±0.20)e—14 

13.1T0.4 

(6.26±0.50)e—14 

15.9T1.3 

(9.19±0.13)e—13 

1.085T0.018 

2007 Jan 24 

(8.17±0.14)e—13 

(6.58±0.18)e—14 

12.4T0.4 

(4.18±0.24)e-14 

19.6T1.2 

(8.04±0.26)e—13 

1.016T0.037 

2010 Feb 27 

(8.89±0.15)e—13 

(7.72±0.68)e—14 

11.5T1.0 

(5.79±0.77)e—14 

15.4T2.1 

(1.64±0.01)e—12 

0.542T0.010 

Table 2. Optical 

emission line fluxes and ratios 






Observation 

[O III] A5007 

[O III] A4959 

[O HI]/ 

[Fe VII] A3759 

[O HI]/ 

[Fe VII] A5159 

[O HI]/ 

Date 

(erg/s/cm 2 ) 

(erg/s/cm 2 ) 

[O III] 

(erg/s/cm 2 ) 

[Fe VII] 

(erg/s/cm 2 ) 

[Fe VII] 


2002 Apr 11 

(8.52±0.03)e—12 

(2.96±0.03)e—12 

2.88T0.03 

(7.87±1.38)e—14 

108.3T19.0 

(4.23±0.57)e—14 

201.4T27.0 

2004 May 28 

(1.18±0.01)e-ll 

(4.60±0.04)e—12 

2.57T0.02 

(1.85±0.25)e—13 

63.8T8.7 

(6.60±0.73)e—14 

178.8T19.7 

2006 Jan 5 

(1.17±0.01)e-ll 

(3.92±0.02)e—12 

2.99T0.02 

(9.66±0.58)e—14 

121.1T7.3 

(5.36±0.21)e—14 

218.3T8.7 

2006 May 29 

(9.80±0.04)e—12 

(3.86±0.03)e—12 

2.54T0.02 

(9.25±1.36)e—14 

106.0T15.6 

(5.34±0.60)e—14 

183.5T20.7 

2007 Feb 9 

(1.46±0.01)e—11 

(5.09±0.04)e—12 

2.87T0.03 

(7.45±1.20)e—14 

196.0T31.6 

(5.32±0.40)e-14 

274.4T20.9 

2010 Feb 18 

(9.79±0.03)e—12 

(3.14±0.02)e—12 

3.12T0.03 

(6.82±1.15)e—14 

143.6T24.3 

(3.46±0.29)e—14 

283.0T24.1 

Observation 

[O III] A5007 

[Fe VII] A5721 

[O III]/ 

[Fe VII] A6087 

[O III]/ 

Ha 

[O III]/ 

Date 

(erg/s/cm 2 ) 

(erg/s/cm 2 ) 

[Fe VII] 

(erg/s/cm 2 ) 

[Fe VII] 

(erg/s/cm 2 ) 

Ha 

2002 Apr 11 

(8.52±0.03)e—12 

(8.67±0.68)e—14 

98.3T7.7 

(1.19±0.10)e—13 

71.6T5.7 

(1.43±0.01)e—11 

0.596T0.005 

2004 May 28 

(1.18±0.01)e-ll 

(9.75±0.81)e—14 

121.0T10.1 

(1.91±0.11)e—13 

61.8T3.6 

(1.47±0.01)e—11 

0.803T0.005 

2006 Jan 5 

(1.17±0.01)e—11 

(8.73±0.39)e—14 

134.0T6.0 

(1.64±0.07)e—13 

71.3T3.2 

(1.64±0.01)e—11 

0.713T0.003 

2006 May 29 

(9.80±0.04)e-12 

(8.32±0.68)e—14 

117.8T9.7 

(1.26±0.11)e—13 

77.8T6.8 

(1.46±0.01)e—11 

0.671T0.005 

2007 Feb 9 

(1.46±0.01)e—11 

(9.93±0.62)e—14 

147.0T9.1 

(1.73±0.08)e—13 

84.4T4.1 

(1.41±0.01)e—11 

1.035T0.007 

2010 Feb 18 

(9.79±0.03)e—12 

(7.25±0.56)e—14 

135.0T10.4 

(1.39±0.09)e—13 

70.4±4.4 

(1.89±0.02)e—11 

0.518T0.005 


following the temporal changes of the coronal lines in rela¬ 
tive rather than absolute flux. In particular, we scale the 
coronal line emission to that of a strong, forbidden low- 
ionisation emission line that is unblended and observed in 
the same spectrum. The emission region that produces the 
low-ionisation narrow lines is believed to be located at large 
enough distances from the central ionising source for its flux 
to remain constant on timescales of decades. We have scaled 
the near-IR and optical coronal lines to the [S III] A9531 and 
[O III] A5007 emission lines, respectively. 

The error estimates on the line fluxes are crucial for as¬ 
sessing the significance of the coronal line variability. The 
data are of relatively high signal-to-noise ratio (continuum 
S/N > 50 — 100) and, therefore, the main sources of mea¬ 
surement errors are of a subjective nature, namely, the place¬ 
ment of the local continuum and related to this the setting 
of the extension of the emission line. These are problematic 
in most cases since the lines in question sit on top of broad 
emission lines, namely, [S III] A9531 in the blue wing of Pae, 
[S VIIl] in the blue wing of Pad, [S ix] in the blue wing of 
Pa/9, [Si Vi] in the red wing of Brd, [O III] A5007 in the red 


wing of H /3 and [Fe x] in the blue wing of Ha. Therefore, in 
order to estimate meaningful uncertainties for the measured 
line fluxes we have considered in addition to the statisti¬ 
cal errors due to the data quality also profile comparisons. 
These were done in velocity space between the scaling lines 
[S III] A9531 and [O III] A5007 and the coronal lines and 
helped us judge the influence of the continuum placement 
on the recovery of the true line profile. The total estimated 
lcr uncertainties for all measured line fluxes are listed in Ta¬ 
bles [T] and [2] They are ~ 1 — 2% for the strongest lines and 
~ 3 — 15% for the weakest ones. 

In order to further constrain the significance of the 
observed coronal line variability, we have considered the 
two extreme cases of where no variability and the high¬ 
est variability are expected. For both scaling lines we ob¬ 
serve also the other emission line that is emitted from the 
same upper level, namely, [S III] A9069 and [O III] A4959. 
Their observed ratios, which should be close to the the¬ 
oretical values of [S ill] A9531 / [S in] A9069= 2.58 and 
[O III] A5007/[0 III] A4959=2.92 (iKramida et al.ll2013lj . are 
not expected to vary and so their observed variability sets a 
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Table 3. Near-IR continuum fluxes and optical coronal line ratios 


Observation 

Date 

A/ l/i,m 

(erg/s/cm 2 ) 

A/l^im / 

[S III] 

A/ 2 .1/xm 

(erg/s/cm 2 ) 

A/2.1 M m/ 

[S III] 

7d U st 

(K) 

[Fe VII] A6087/ 
[Fe VII] A3759 

[Fe VII] A5159/ 
[Fe VII] A6087 

2002 Apr 23 

(8.57±0.03)e—11 

112.6T1.5 

(1.22±0.01)e—10 

160.3T1.9 

1316 

1.512T0.290 

0.355T0.055 

2004 May 23 

(7.79±0.25)e—11 

73.5T2.9 

(1.18±0.01)e—10 

111.3±2.7 

1281 

1.032T0.153 

0.346T0.043 

2006 Jan 8 

(9.34±0.07)e—11 

79.8T0.7 

(1.62±0.01)e—10 

138.5T1.0 

1328 

1.698T0.127 

0.327T0.020 

2006 Jun 12 

(6.67±0.07)e—11 

66.9T1.0 

(1.50±0.01)e—10 

150.5T1.7 

1278 

1.362T0.233 

0.424T0.061 

2007 Jan 24 

(5.12±0.28)e—11 

62.7T3.6 

(1.04±0.01)e—10 

127.3T2.4 

1231 

2.322T0.390 

0.308T0.028 

2010 Feb 27 

(2.18±0.01)e—10 

245.2T4.4 

(3.47±0.01)e—10 

390.3T6.6 

1393 

2.038T0.367 

0.249T0.026 


Table 4. XMM-iVerofon archival observations 


No. 

Observation 
date and time 

Exp. 

(ks) 

Obs ID 

Campaign 

1 

2000:12:21 16:45:30 

33 

0112310101 

2000 Dec 

2* 

2000:12:21 14:21:13 

8 

0112310501 


3 

2000:12:22 10:42:18 

62 

0112830201 


4 

2000:12:22 02:53:42 

23 

0112830501 


5* 

2000:12:22 04:02:18 

5 

0112830601 


6 

2003:05:25 08:37:56 

19 

0143500101 

2003 May 

7 

2003:05:26 20:35:02 

19 

0143500201 


8 

2003:05:27 15:16:51 

19 

0143500301 


9 

2006:05:16 06:21:42 

40 

0402660101 

2006 May 

10 

2006:11:29 17:20:13 

53 

0402660201 

2006 Nov 

11 

2011:05:11 15:20:09 

16 

0657840101 

2011 May 

12 

2011:06:12 13:36:34 

16 

0657840201 



* no EPIC data 


lower threshold for the significance of the coronal line vari¬ 
ability. Then, we have measured the fluxes of the two promi¬ 
nent broad emission lines Pa /3 (in the near-IR) and Ha (in 
the optical). Since the BELR is expected to be the most 
variable of any AGN emission line region, the observed vari¬ 
ability of these broad lines gives an estimate of the maximum 
value that can be reached within the current data set. These 
results are also listed in Tables □ and [2] 

In addition to the emission lines, we have measured in 
the near-IR spectra the continuum fluxes at the rest-frame 
wavelengths of ~ 1 /im and ~ 2.1 ) im (see Table [3j). As 
we have shown in lLandt et al.l (l201ll ) , the former is dom¬ 
inated by the accretion disc flux, which is believed to be 
the main source of ionising radiation in AGN and so the 
driver of the observed variability, whereas the latter is emit¬ 
ted from the hot dust component of the obscuring torus, 
which, if its location indeed coincides with that of the coro¬ 
nal line region, should have a variability response similar 
to it. Furthermore, we have derived the hot dust tempera¬ 
ture from blackbod y fits to the near-IR spectral continuum 
as described in lLandt et ali rt201ll ') and list it also in Table 
□ In Section □ we will compare these values to the op¬ 
tical coronal line ratios [Fe VIl] A6087/[Fe VIl] A3759 and 
[Fe VI 1 ] A5159/[Fe vn] A6087, whic h are suitable indicators 
of the gas temperature and d ensity jNussbaumer et alii 19821 ; 
[Keenan fc Norringtonlll 9871 1. respectively (listed in Table[3]). 


2.2 X-ray spectroscopy 

Between December 2000 and June 2011, the time period that 
overlaps with our optical and near-IR observations, there 


were five observational campaigns with XMM-Newton on 
NGC 4151 (see Table [p. We have reduced these archival 
data sets using the standard tasks in the XMM-Newton 
SAS software (version 13.5.). All the X-ray spectral anal¬ 
ysis was done using the SPEX software^ dKaastra et ali 
1996) . We have a ssume d solar abundances as given by 
Lodders fc Palmel d2009l) and corrected all spectra for a 
Galactic column density of 2xl0 24 m 2 . Three of the cam¬ 
paigns have multiple observations, with two of them (2000 
December and 2003 May) extending over a period of 2-3 
days and another one (2011 May) spanning a period of a 
month. Since the X-ray emissio n lines are not expected to 
vary on days to weeks timescales dPetmers et al.ll2008l . 120091) 
and in order to improve the S/N ratio, we have stacked the 
individual RGS spectra into one spectrum per campaign. 
We have treated the May 16 and November 29, 2006, ob¬ 
servations separately, since the time span between them is 
about half a year. The longest and shortest observations are 
for the December 2000 and May 2011 campaigns with total 
exposure times of 131 ks and 32 ks, respectively. 

In the X-rays, NGC 4151 is classified as an obscured 
Seyfert 1 galaxy, i.e. the soft X-ray emission is heavily ob¬ 
scured by an absorber. However, this absorption is not as 
deep as that observed for Seyfert 2 galaxies and a very 
weak continuum is present. This absorber is different from 
the warm absorber detected in the ultravio let (UV) and X- 
ray s pectra of ~50% of Seyfert 1 galaxies (ICrenshaw et al.1 
2003). It has a much larger column density, outflow veloc¬ 
ity and velocity broadening, is variable on much shorter 
timescales and only partially covers the X-ray continuum 
source. In order to d ifferentiate t hi s abso rber from the warm 
absorber, we follow [Kaastra et al.l (120141 s ) and will call it ’the 
obscurer’. Since the soft and hard X-ray continuua as well 
as the properties of the obscurer, such as its column den¬ 
sity, covering factor and potentially also the ionisation, are 
known to be variable on day timescales or less, we have not 
combined the EPIC pn spectra, but fitted them separately. 

The presence of the obscurer in NGC 4151 allows us 
to study in detail the X-ray emission lines, which would 
otherwise be swamped by the high continuum flux. How¬ 
ever, the obscurer also complicates the analysis, in particular 
the continuum fitting, since the absorption is so deep that 
the individual absorption lines blend and create a pseudo¬ 
continuum. Decomposing the total absorption into obscurer 
and warm absorber and characterising each individually is 
not possible without prior knowledge of the ionisation struc¬ 
ture of the warm absorber. This, however, is not known for 

1 http://www.sron.nl/spex 


© 2015 RAS, MNRAS 000.HHT21 






































Coronal line variability in NGC 4151 5 


Table 5. X-ray emission line and continuum fluxes 


Observation 

Epoch 

O VII f 0.561 keV 
(ph/s/cm 2 ) 

O VII i 0.569 keV 
(ph/s/cm 2 ) 

O VII r 0.574 keV 
(ph/s/cm 2 ) 

O VIII a 0.654 keV 
(ph/s/cm 2 ) 

Ne IX f 0.905 keV 
(ph/s/cm 2 ) 

2000 Dec 

2003 May 
2006 May 
2006 Nov 

2011 May 

(5.05±0.17)e—04 
(5.10±0.21)e—04 
(4.704=0.21)e-04 
(3.89±0.24)e—04 
(5.63=t0.95)e—04 

(0.70±0.09)e—04 
(0.434=0.09)e-04 
(0.784=0.ll)e—04 
(0.554=0.31)e—04 
(0.984=0.58)e—04 

(1.514=0.10)e—04 
(1.364=0.13)e—04 
(1.554=0.13)e—04 
(1.494=0.34)e-04 
(1.804=0.62)e—04 

(1.474=0.06)e—04 
(2.02=t0.06)e—04 
(1.594=0.08)e-04 
(1.604=0.08)e-04 
(1.874=0.28)e-04 

(0.594=0.04)e—04 
(0.564=0.04)e-04 
(0.454=0.08)e-04 
(0.494=0.07)e-04 
(0.874=0.16)e-04 

Observation 

Epoch 

N VI f 0.420 keV 
(ph/s/cm 2 ) 

N VII a 0.500 keV 
(ph/s/cm 2 ) 

C VI a 0.368 keV 
(ph/s/cm 2 ) 

/2-10keV 

(erg/s/cm 2 ) 


2000 Dec 

2003 May 
2006 May 
2006 Nov 

2011 May 

(1.51±0.09)e—04 
(1.594=0.10)e—04 
(1.274=0.13)e—04 
(1.234=0.14)e-04 
(1.054=0.48)e-04 

(0.87±0.05)e—04 
(0.974=0.06)e—04 
(0.94±0.08)e—04 
(0.734=0.09)e—04 
(0.844=0.28)e-04 

(2.664=0.13)e-04 
(2.914=0.15)6-04 
(3.104=0.22)e—04 
(2.434=0.22)e—04 
(2.634=0.74)e—04 

(6.374=0.06)e—11 
(3.05=t0.24)e—10 
(7.914=0.08)e-ll 
(1.374=0.01)e-10 
(1.554=0.05)e-10 



NGC 4151 since it has never been observed in an unob¬ 
scured state. However, for the Seyfert 1 galaxy NGC 5548 
we have high-quality X-ray spectra in both the un obscured 
llSteenbrueee et aDl2003l. 2005^ Detmers et al.ll200^ 1 and ob¬ 
scured state ( Kaastra et al.ll2014T ). Thus, we have assumed 
in our X-ray analysis that the obscurer in NGC 4151 has the 
same structure as that in NGC 5548, albeit being deeper. 
This is supported by the fact that the C IV absorption 
lines observed in the Hubble Space Teles cope (HST) UV 
spectra of NGC 4151 llKraemer et al.ll20()d ) and NGC 5548 
([Kaastra et al.li2014l ) have a similar width and blueshift. 

Following the approach chosen bv lKaastra et all ( 2014 1 
for NGC 5548, we first fitted the May 2003 spectrum, since 
it is the least absorbed state and allows for a better deter¬ 
mination of the warm absorber and obscurer components. 
Each warm absorber or obscurer component was fitted with 
an xabs model. Then, using their results for NGC 5548, we 
assumed that the obscurer has two components, one close to 
or neutral with a lower covering factor and another medium- 
ionised with a much larger covering factor. From the UV 
spectrum of NGC 4151 we know that there are at least six 
war m absorber com ponents with different velocities (see Fig. 
2 in lKraemer et al.ll2006l l. but their ionisation structure is 
undetermined by these data. We have added a further four 
warm absorber components (each an xabs model) with dif¬ 
ferent ionisation parameters. Although the ionisation pa¬ 
rameter of the warm absorber and obscurer can overlap, 
we only allowed higher ionization parameters for the warm 
absorber. The free parameters in the fit were the hydrogen 
column density and ionisation parameter, and in the case 
of the obscurer also the covering factor. The outflow veloc¬ 
ity and velocity broadening were left at the standard values 
since they cannot be constrained. Finally, we assumed that 
the narrow emission lines are unaffected by the obscurer and 
warm absorber. 


We modelled the continuum as observed with EPIC-pn 
with a reflection component and a comptoionisation model, 
which dominate the hard X-rays, and added a modified 
blackbody in the soft X-ray regime. The fit also includes, 
but the parameters are frozen, the best fit model for the 
emission lines, both narrow and broad, and the radiative re¬ 
combination continua as determined from the RGS spectra. 
The Fe Ka and K/3 lines, which can only be studied with 


the EPIC instruments, were fitted as Gaussians. Both lines 
are unresolved. We fit the spectra in the 0.3— 10 keV energy 
range and used the option of optimal binning available in 
SPEX. This gives an acceptable fit to the EPIC-pn spectra. 
However, due to the heavy absorption, there is a certain de¬ 
generacy between the normalisation and temperature of the 
modified blackbody, the hydrogen column density, ionisa¬ 
tion parameter and for the obscurer also the covering factor, 
of the absorption components and the normalisation of the 
emission lines and radiative recombination continua (frozen 
in the EPIC pn fit). In Tabled we list the fluxes and lcr er¬ 
rors for the eight strongest X-ray coronal lines, mainly from 
oxygen but also from nitrogen, carbon and neon, and the un¬ 
absorbed 2 — 10 keV continuum fluxes. The errors on the line 
fluxes are ~ 3 — 10% for the strongest lines and ~ 10 — 20% 
for the weakest lines. The errors for the May 2011 observing 
campaign are much larger on all measurements (~ 20—60%). 


3 THE VARIABILITY BEHAVIOUR 

3.1 The near-IR and optical coronal lines 

In the ~ 8 year period sampled by the near-IR and optical 
spectroscopy the ionising flux of the accretion disc, which 
is assumed to be the main driver for the variability of the 
broad line, coronal line and hot dust emission regions, de¬ 
creased by a factor of ~ 2 in the first five years and increased 
by a factor of ~ 4 in the following three years (see Fig. [I]). 
We note that the flux increase of ~ 9% observed between 
the May 2004 and January 2006 epochs is consistent with 
zero once the deviations from the theoretical value of the 
[S ill] A9531/[S ill] A9069 ratios measured in these spectra 
are taken into account (~ 5.5% and ~ 2.9%, respectively). 
The Pa/1 broad line region and hot dust emission responded 
to the accretion disc variability in a similar way (see Figs. 
[2]and[T] respectively); their flux decreased in the first two 
years by ~ 70% and ~ 50%, respectively, and then continu¬ 
ously increased for the next six years by a factor of ~ 3 and 
~ 4, respectively. We note that the flux decrease of ~ 20% 
observed for the hot dust emission between the June 2006 
and January 2007 epochs is similar to the deviation from the 
theoretical value measured for the [S III] A9531/ [S III] A9069 
ratio in the latter spectrum (~ 13.3%). The flux of the Ha 
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Figure 1. The variability of the continuum fluxes at rest-frame 
vavelengths of ~ 1 /im (sampling the accretion disc), ~ 2.1 /im 
sampling the hot dust emission) and in the energy range 2 — 
L0 keV (related to the accretion disc). The near-IR continuum 
luxes were divided by the [S III] A9531 line flux and all continuum 
luxes were normalised to the value of the May/June 2006 epoch. 
vVe plot lcr error bars. 
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Figure 2. The variability of the near-IR coronal lines in the pe¬ 
riod April 2002 - February 2010. For comparison, we show also 
the [S III] A9069 narrow line and Pa/3 broad line, which are 
expected not to vary and to vary maximally, respectively. The 
[S III] A9531/[S III] A9069 flux ratios were normalised to the the¬ 
oretical value, whereas those between [S III] A9531 and the other 
lines were normalised to the value of the June 2006 epoch. We 
plot lcr error bars. The ionisation potentials of the coronal lines 
are given at the top right. 


broad line region also decreased in the first two years (by 
~ 35%) and increased in the next two years (by ~ 20%), 
but then it decreased for a short period (between June 2006 
and January 2007 by ~ 54%) before it increased strongly 
for the last three years (by a factor of ~ 2; see Fig. 0. The 
intermittent flux decrease might be caused by a reddening 
event, which we discuss further in Section 0] 

The variability response of the near-IR coronal lines in 
the four year period between January 2006 and February 
2010 is similar to that of the Pa/3 broad line region and hot 
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Figure 3. Same as Fig. [2] for the optical coronal lines. For com¬ 
parison, we show also the [O III] A4959 narrow line and Ha broad 
line, which are expected not to vary and to vary maximally, re¬ 
spectively. 
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Figure 4. The variability of the X-ray coronal lines in the period 
December 2000 - May 2011. The flux values were normalised to 
the value of the May 2006 epoch. We plot lcr error bars. The 
ionisation potentials of the coronal lines are given at the top right. 


dust emission, with a trend that the higher the ionisation 
potential, the higher the flux increase, namely, ~ 69% for 
the [Si x] line, ~ 28% for the [S ix] line, ~ 38% for the 
[S Vlll] line, and no significant variability for the [Si Vi] line 
(see Fig. n. In the first two years sampled by the data sig¬ 
nificant variability is observed only for the [S VIIl] line; a 
flux increase by ~ 54%. Then, for all but the [S ix] line, a 
flux decrease (instead of the increase seen for the broad line 
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region and hot dust emission) is observed between May 2004 
and January 2006, namely, ~ 52% for the [Si x] line, ~ 55% 
for the [S VIIl] line, and ~ 28% for the [Si Vi] line. 

Between April 2002 and January 2006 we detect signif¬ 
icant variability for only one optical coronal line, namely, 
[Fe VIl] A3759, which behaves similar to the [S VIIl] line (see 
Fig. [3]); a flux increase (by ~ 70%) followed by a flux de¬ 
crease (by ~ 90%). The flux decrease observed for the Ha 
broad line region between June 2006 and January 2007 is 
clearly evident for the [Fe VIl] A3759, [Fe VIl] A5159 and 
[Fe VIl] A5721 lines; their fluxes changed by ~ 85%, ~ 50% 
and ~ 25%, respectively. If we consider the period of Jan¬ 
uary 2006 to January 2007 instead, we observe a similar 
behaviour also for the [Fe VIl] A6087 line; its flux decreased 
by ~ 18%. In the following three years the [Fe VIl] A6087 line 
flux increased by ~ 20%, as observed for the broad line re¬ 
gion and hot dust emission, whereas none of the other optical 
coronal lines varied signihcantly. 


3.2 The X-ray coronal lines 

In the ~ 10.5 year period sampled by the X-ray spectroscopy 
the unabsorbed 2—10 keV continuum flux, which is assumed 
to be produced by the central ionising source and so to be 
linked to the accretion disc flux, increased by a factor of ~ 5 
in the first two and a half years, decreased by a factor of 
~ 4 in the next three years, and then increased again by a 
factor of ~ 2 in the next five years. 

The variability response of the three forbidden X-ray 
coronal lines in the period between November 2006 and May 
2011 is similar to that of the near-IR coronal lines and that 
of the Pa/3 broad line region and hot dust emission, and 
again with a trend that the higher the ionisation potential, 
the higher the flux increase, namely, a factor of ~ 2 for the 
Ne IX f line, ~ 45% for the O VII f line and no significant 
variability for the N VI f line (see Fig. [4]). Then, whereas the 
flux stayed roughly constant in the first two and half years 
sampled by the data, similar to the near-IR coronal lines, 
a flux decrease (instead of the increase seen for the broad 
line region and hot dust emission) is observed in the period 
between May 2003 and May 2006 for the O VII f line (by 
~ 31%) and the N VI f line (by ~ 30%). 

In the entire time period sampled by the data no sig¬ 
nificant variability is observed for the recombination line 
O VII r, whereas the intercombination line O VII i showed 
strong variability between December 2000 and May 2006, 
with a flux decrease by ~ 62% followed by a flux increase 
by ~ 82%. No significant variability is observed for the 
O VII i line after May 2006, which could be due to the rel¬ 
atively large errors on the line flux for the following two 
observing epochs. The O VIII a line showed also significant 
variability between December 2000 and May 2006, but con¬ 
trary to the O VII i line, its flux increased first (by ~ 37%) 
and then decreased (by ~ 27%). Also for this line there 
is no significant variability after May 2006. The variabil¬ 
ity response of the two permitted X-ray coronal lines with 
the highest ionisation potentials is similar to that of the 
O VII f line; a significant flux decrease is observed for the 
N VII a line only between May 2003 and November 2006 (by 
~ 34%) and for the C VI a line only between May 2006 and 
November 2006 (by ~ 28%). 


4 THE ORIGIN OF THE CORONAL LINE 

EMISSION REGION 

The current understanding is that the coronal line region 
in AGN is dust-free, photoionised and located beyond the 
BELR at distances from the central ionising source similar 
to those of the hot inner face of the obscuring dusty torus. 
In this scenario, we expect to observe the coronal line flux 
to vary similarly to that of the broad lines and hot dust 
emission but on timescales longer and shorter than the for¬ 
mer and latter, respectively. With our data set we can test 
the first premise, although its sparse time sampling does not 
allow us to measure any variability timescales. 

The variability behaviour observed for the broad lines 
and hot dust emission in the period covered by our data is 
both similar and simple; a flux decrease of ~ 40 — 70% in 
the first two years followed by a strong flux increase by a 
factor of ~ 2 — 4 in the following six years. This variability 
behaviour is not observed for any of the coronal lines. In 
general, the coronal lines varied weakly if at all, with the 
largest flux change observed to be only ~ 50 — 90% in a 
period of about two years. Specifically, in the first four years 
sampled by the data, the coronal lines either did not vary 
significantly or showed the opposite behaviour to that of 
the broad lines and hot dust emission. In the last four years 
sampled by the data, only the coronal lines with the highest 
ionisation potentials showed a variability behaviour similar 
to that of the broad lines and hot dust emission, but with a 
much reduced amplitude (a lower flux change per year by a 
factor of ~ 2 — 4), whereas the flux of the coronal lines with 
relatively low ionisation potentials remained unchanged. 

The characteristic variability response time is a sum of 
the light travel time, which depends mainly on the location 
of the emission region, and the recombination time, which 
depends strongly on the gas number density: 

Tvar — Tit T Tree (1) 

Therefore, the low variability amplitude observed for the 
coronal line gas in NGC 4151 indicates either that it is lo¬ 
cated well beyond the broad line region and dusty torus or 
that it has a relatively low density or both. In the follow¬ 
ing, we use published results from high-spatial resolution 
imaging campaigns at near-IR, optical and X-ray frequen¬ 
cies and apply plasma diagnostics to our own data in order 
to constrain both the location and gas number density of 
the coronal line emission region. 


4.1 Published high-spatial resolution observations 


NGC 4151 has a well-me asured h y droge n broad line region 
lag time of about a week |]Zn et ahll201ll') and a hot dust la g 
time varying between ~ 30 — 70 days ( Koshida et al. I [20141 ). 
One possible explanation for the fact that the coronal line 
region reacted much weaker to changes in the ionising flux 
than both the broad line and hot dust emission regions is 
that it is located much further out than them. 

High-spatial resolution observations of NGC 4151 
in _the_ near-IR were prese nted by several authors. 
IStorchi-Bergmann et al.l (|2009| J used the Gemini Near-IR 
Integral Field Spectrograph (NIFS) on December 2006 to 
image the source in the 0.94 — 2.42 /am wavelength range 
with a spatial resolution of 0.12 — 0.16 arcsec, which corre- 
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sponds to « 8 — 10 pc at the source. Their contour maps of 
the [S VI11] and [S ix] coronal li nes s ho w a compact , spa¬ 
tially unresolved reg ion. iMuller-Sanchez et al.l (120111 1 and 
llserlohe et al.l (120131 1 used the Keck OH Suppressing In- 
fraRed Imaging Spectrograph (OSIRIS) on March 2006 and 
February/May 2005, respectively, to image the source in the 
1.96 — 2.38 p m wavelength range with a spatial resolution 
of 0.08 — 0.11 arcsec, which corresponds to ~ 5 — 7 pc at 
the source. Their contour maps of the [Si Vi] coronal line 
show a bright nucleus containing ~ 70% of the total emis¬ 
sion, which is surrounded by low-level extended emission up 
to « 50 — 80 pc. 

The source NGC 4151 was observed with the Space 
Telescope Imaging Spectrograph (STIS) on-board the Hub¬ 
ble Space Telescope (HST) on several occasions. These ob¬ 
servations span a large optical wavelength range and have a 
spatial resol ution of 0.1 arcsec , corresponding to « 7 pc at 
the source. iNelson et al.l (|2000l l presented measurements of 
the [Fe vn] A3759, [Fe vn] A5721 and [Fe vn] A6087 coronal 
line fluxes along the slit at two different position angles for 
observations taken between January and June 1998. Along 
one of the position angles (P.A. 70°), the [Fe VIl] A3759 emis¬ 
sion was spatially unresolved, but the [Fe VIl] A5721 and 
[Fe vnj A6087 emission lines showed extents of up to 
sa 70 — 100 pc from the nucleus. Along the other position an¬ 
gle (P.A. 221°), all three [Fe VIl] emission lines were spatially 
resolved, however, with much lesser extents of « 20 — 50 pc 
from the nucleus. In all cases, the total [Fe vn] emission was 
do minated b^_the nuclea r flux. 

[Wang etahl (l2011bl lah analysed the deep exposure of 
NGC 4151 taken in March 2008 with the Advanced CCD 
Imaging Spectrometer (ACIS) on-board the Chandra X-ray 
observatory. By applying subpixel event repositioning and 
binning techniques they were able to improve the effective 
spatial resolution of the ACIS images to better than 0.4 arc¬ 
sec, which corresponds to ss 25 pc at the source. Their emis¬ 
sion line maps of O VII (f, i and r), O VIII a and Ne IX (f, i 
and r) show an extremely bright nucleus surrounded by ex¬ 
tended emission in the inner region of « 130 pc, and for the 
O VII and O VIII emission lines also very low-level emission 
(~ 10% of the total flux) up to « 2 kpc. 

In summary, high-spatial resolution observations of 
NGC 4151 show that whereas its coronal line emission re¬ 
gion is extended in some chemical species, the total flux that 
our variability study is most sensitive to is dominated by the 
unresolved nuclear region. This region can currently be con¬ 
strained by direct imaging only to < 5 pc (or < 15 light 
years), which is much larger than the measured distance of 
the hot dust emission of ~ 2 light months. 

4.2 Plasma diagnostics 

For the optical and X-ray coronal line gas we can esti¬ 
mate the density using the three iron lines [Fe VIl] A3759, 
[Fe VIl] A5159 and [Fe VIl] A6087 and the three lines 
from the helium-like ion of oxygen O VII f, O VII i and 
O VII r, respectively. Usin g the plasma simul ation code 
CLOUDY (last described by Ferland et al.l 1199a) . we have 
generated a temperature versus density grid for the col- 
lisionally excited [Fe VIl] lines (see Fig. 0 and a ioni¬ 
sation parameter versus density grid for the O VII lines 
(see Fig. 0. The iron collision strengths were taken from 



Figure 5. The observed optical coronal line ratios 
[Fe VII] A6087/A3759 versus [Fe VII] A5159/A6087, overlaid 
with curves of constant temperature (logT = 4,4.1,4.2, ...5 K) 
and constant number density (logn e = 3, 3.2, 3.4, ...7 cm 3 ) for 
the case of photoionisation equilibrium. The case of collisional 
ionisation equilibrium is shown in green. We plot I a error bars. 


IWitthoeft fc Badnell 1200 81. For the photoionisation simu¬ 
lations we have approximated the incident radiation field 
with th e mean AGN spectral en ergy distribution de¬ 
rived bv [Mathews fc Ferlandl (1987). As discussed by pre¬ 
vious studies, the line ratios [Fe VIl] A6087/A3759 and 
[Fe vn] A5159/A6087 are suitab le indicators of temper¬ 
ature and density, respect ively (Nussbau mer et all Il982l : 
iKeenan fe Norringtonlll987l ). whereas the X-ray line ratios 
O VII (f+i)/r (the so-called G ratio) and O VII f/i (the so- 
called R ratio) trace the ionisation parameter, which is di¬ 
rectly related to the kinetic gas temper ature, and densi ty, 
respectively, for a given colum n density dPorauet fe Dubaul 
120001 : iPorter fe Ferlandl l2007fl . We assumed the cases of ei¬ 
ther photoionisation or collisional ionisation equilibrium. 

The measurements for NGC 4151 in both Figs. [5] and 
0 give two main consistent results; the coronal line gas is 
photoionised rather than collisionally ionised and its density 
appears to be relatively low. In the case of the optical coronal 
line gas, if we assume that the gas density did not change 
with time, four of the observing epochs constrain it to n e ~ 
10 3 cm -3 within ~ lcr, whereas the two observing epochs 
January and June 2006 reach this value within ~ 2a. If we 
assume the gas density changed between observing epochs, 
at the 2 ct level it varied between n e ~ 10 5 ’ 2 cm~ 3 and n e < 
10 3 cm -3 , since at the low-density limit the grid contours 
fall on top of each other. The highest density value is about 
two orders of magnitude lower than the critical density of 
[Fe VI 1 ] of ricrit ~ 3 x 10 7 cm -3 . We get a similar result 
for the X-ray coronal line gas; assuming the gas density did 
not change with time, two observing epochs (November 2006 
and May 2011) constrain it to n e ~ 10 3 cm -3 within ~ la, 
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Figure 6. The observed X-ray coronal line ratios O VII (f+i)/r (G 
ratio) versus O VII f/i (R ratio), overlaid with curves of constant 
ionisation parameter (log U = — 1, — 0.5, 0, ...2) for two constant 
number densities (n e = 10 3 cm -3 and n e = 10 12 cm -3 ) and a 
fixed column density of Nn = 10 21 cm -2 , assuming photoionisa¬ 
tion equilibrium. The case of collisional ionisation equilibrium is 
shown in green. We plot 1 a error bars. 


whereas three observing epochs reach this value within ~ 3a. 
If we assume the density changed between observing epochs, 
at the 2a level it varied between n e ~ 10 12 cm~ 3 (i.e. close 
to the critical density) and n e < 10 3 cm' 3 , since at the 
low-density limit the grid contours fall on top of each other. 

It is worth noting that all measured O VII ratios and 
most of the [Fe VIl] ratios are to the right of the theoretical 
grids. In the case of O VII, there are no known uncertainties 
in the atomic data and adding a contribution from collision- 
ally ionised plasma would not increase the R ratios, which 
are observed to be much higher than predicted by theory. 
Similarly high R rati os are found also in o ther AGN, e.g., 
R = 5.50 for Mrk 3 llP ounds fe Paeel [20051 ) and R = 7.73 
for NGC 1068 llKraemer et al.ll2015l ). However, in the case 
of [Fe VIi], it is possible that some contribution from flu¬ 
orescence is present given the relatively low gas density in 
combination with the large ionisation parameter we esti¬ 
mate below. The large number of permitted Fe VII lines 
in the far-UV can be pumped by the continuum, which, 
when radiatively decaying, will add to the lower levels re¬ 
sponsible for the [Fe VIl] lines studied here. Fig. [7] shows 
the importance of this effect for a gas of number density 
n e = 10 3 cm -3 in dependence of the ionisation parameter. 
At a ionisation parameter of log U ~ 1 the expected increase 
in the [Fe vn] A5159/A6087 ratio is a factor of ~ 1.5, whereas 
the increase in the [Fe VIl] A6087/A3759 ratio is only ~ 10%. 

From Fig. 0 we estimate a temperature of T ~ 
10 4 -25 18000 K for the Fe 6+ gas. At this temperature, 

the recombination time is: 


Tree 


1 

He ' 1.70 X 10- 10 S ’ 


( 2 ) 



log U 


Figure 7. The effect of continuum pumping on the optical 
coronal line ratios [Fe VII] A6087/A3759 (red solid line) and 
[Fe VIl] A5159/A6087 (black solid line) in dependence of the ion¬ 
isation parameter. The gas is assumed to have a temperature of 
T e = 10 4-25 K and a number density of n e = 10 3 cm -3 . 



Figure 8. The density run with radius calculated for the rela¬ 
tionship n oc r ~ 2 (black solid line) and scaled to the observed lag 
time and expected gas number density for the hydrogen broad line 
region in NGC 4151 (black dashed lines). The red dashed lines 
show the observed lag time and expected gas number density for 
the hot dust. The estimated density and distance of the coronal 
line region are indicated by the blue dashed lines. 


which for a gas number density of n e = 10 3 cm -3 results 
in Tree ~ 5.9 x 10 6 s ~ 2.2 months. We note that we 
have used the total recombination coefficient, i.e. the sum 
of the radiative and dielectronic recombination coefficients, 
wh ich was c a lculat ed based on the atomic data presented 
by iGu et all (|2006l ). This value differs from that listed in 
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Figure 9. The top and bottom panels show the continuum flux 
at rest-frame wavelength of ~ 1 /im (sampling the accretion disc) 
and the optical coronal line ratio [Fe VII] A6087/A3759, which is 
a suitable indicator of the gas temperature (the higher its value, 
the lower the temperature; see Fig. 0 , respectively, versus the 
hot dust blackbody temperature. 

lOsterbrock fe Ferlandl (|2006l f by about an order of magni¬ 
tude, since the latter represented only the radiative recombi¬ 
nation coefficient. From Fig. 0 we constrain the ionisation 
parameter for all X-ray observing epochs at the la level 
to log U ~ 1. Then, using the unabsorbed X-ray luminos¬ 
ity of Z/ 2 -iokev ~ 10 43 erg s~ 4 resulting from our fits for 
the highest-flux epoch as a proxy for the ionising luminos¬ 
ity producing O VII and the best-fit X-ray spectral slope of 
this epoch of F ~ 1.43 to estimate the mean ionising pho¬ 
ton energy, the calculated distance of the coronal line region 
from the central ionising source for the above gas density 
is rit ~ 730 light days ~ 2.0 light years. This value, which 
is far below the resolution provided by current direct imag¬ 
ing results (see Section EE), puts this region well beyond 
the hot inner face of the obscuring dusty torus. Together 
our estimates of r rec and rit give a characteristic variabil¬ 
ity response time of r var « 2.2 years, which is consistent 
with our main observational result that the coronal line gas 
emission in NGC 4151 has varied in an observing period of 
~ 8 — 11 years, but only weakly so relative to the variability 
observed for the broad lines and hot dust emission. 

But how can we explain our other observational result 
that during half the time period covered by the data the 
variable coronal lines showed the opposite trend to the flux 
changes observed for the broad lines and hot dust emission, 
i.e. their emission decreased instead of increased? In order 
to answer this question, we first need to understand how the 
coronal line region in NGC 4151 relates to the BELR, the 
dusty torus and the low-ionisation NELR. For this purpose, 
we have plotted in Fig.[H]the run of gas number density with 
radius for a relationship n oc r~ 2 , which keeps the ionisa¬ 
tion parameter constant. We have scaled this relationship 


to the observed lag time for the hydrogen BELR (of about 
a week) and the typical density for this emission region of 
n e ~ 10 10 cm -3 . This scaling gives at the observed lag time 
for the hot dust (of ~ 50 days) a density of n e ~ 10 8 cm~ 3 , 
which is consistent with what is currently assumed for the 
dusty torus material. From the relationship in Fig. [S] we 
see that, at the density estimated for the coronal line re¬ 
gion of n e ~ 10 3 cm -3 , a ionisation parameter typical of the 
BELR is reached at distances from the ionising source of 
~ 55 light years, which is in the range observed for the spa¬ 
tially resolved low-ionisation NELR. However, the coronal 
line region in NGC 4151 appears to require a much higher 
ionisation parameter than the hydrogen BELR for its lines 
to form, which, at a given radius, can only be achieved by 
a considerable drop in density. Therefore, it seems that the 
coronal line region is not part of a continuous gas distribu¬ 
tion but rather an independent entity. 

One possibility proposed by lPier fe Voitl (Il995l i is that 
the coronal line region is a layer on the inner part of the 
dusty torus that becomes an efficient coronal line emitter 
only when evaporated in an X-ray heated wind. If this sce¬ 
nario applies to NGC 4151, the wind will have undergone 
adiabatic expansion from its launch location at the inner 
face of the torus until the gas density and distance from 
the ionising source are optimal to give the required high 
ionisation parameter. In the process, the coronal line gas 
will have cooled. In Fig. 0 we present tentative evidence 
for this scenario. In the bottom panel, we compare the tem¬ 
perature of the coronal line gas as measured by the line 
ratio [Fe VIl] A6087/[Fe VIl] A3759 (the higher its value, the 
lower the gas temperature; see Fig. EJ with that of the hot 
dust. We find that the two temperatures behave in opposite 
ways; the temperature of the coronal line gas is high when 
the hot dust temperature is low and vice versa. Only the 
data from January 2007 is an exception to this trend and 
might indicate a reddening event, which could also explain 
the decreased flux of the Ha broad line relative to that of 
the Pa/3 broad line for this epoch (see Section 13.111 . In the 
top panel of Fig. 0 we show the continuum flux at rest- 
frame wavelength of ~ 1 fi m, which samples the accretion 
disc luminosity, versus the hot dust temperature. A clear 
correlation is apparent, which indicates that the change in 
temperature for the hot dust is due to direct heating by the 
central ionising source. Therefore, the increased AGN ra¬ 
diation that heats the dusty torus appears to increase the 
coo li ng effic i ency o f the coronal line gas. In the scenario of 
Pier fe Voit (Il995f ). the dusty clouds will be evaporated in 
an X-ray heated wind more efficiently for higher AGN lumi¬ 
nosities, which will lead to an increase in mass outflow rate 
but also to a stronger adiabatic expansion and so cooling. 


5 SUMMARY AND CONCLUSIONS 

We have presented the first extensive study of the coronal 
line variability in an AGN. Our data set for the nearby, 
well-known source NGC 4151 is unprecedented in that it 
includes six epochs of quasi-simultaneous optical and near- 
IR spectroscopy spanning a period of ~ 8 years and five 
epochs of X-ray spectroscopy overlapping in time with it. 
Our main results are as follows. 

(i) The variability behaviour observed for the broad 
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emission lines and hot dust emission was not mirrored in any 
of the coronal lines. The coronal lines varied only weakly, if 
at all. Specifically, in the first four years sampled by the 
data, the coronal lines either did not vary significantly or 
showed the opposite behaviour to that of the broad lines 
and hot dust emission, whereas after that only the coronal 
lines with the highest ionisation potentials showed a vari¬ 
ability behaviour similar to that of the broad lines and hot 
dust emission, but with a much reduced amplitude (a lower 
flux change per year by a factor of ~ 2 — 4). 

(ii) We have applied plasma diagnostics to the optical 
[Fe VIl] and X-ray O VII emission lines in order to constrain 
the gas number density, temperature and ionisation param¬ 
eter of the coronal line region. We find that this gas has 
a relatively low density of n e ~ 10 3 cm -3 and requires a 
relatively high ionisation parameter of log U ~ 1. 

(iii) We estimate the distance of the coronal line region 
in NGC 4151 from the central ionising source for the above 
gas density to be rit ~ 2.0 light years. This value, which is 
well below the spatial resolution provided by current direct 
imaging results, puts this region well beyond the hot inner 
face of the obscuring dusty torus (of ~ 2 light months). 
Together with the recombination time this results in a char¬ 
acteristic variability time scale of r va r ~ 2.2 years, which is 
consistent with our main observational result that the coro¬ 
nal line gas emission has varied in an observing period of 
~ 8 — 11 years, but only weakly so relative to the variability 
observed for the broad emission lines and hot dust emission. 

(iv) Since the coronal line region requires a much higher 
ionisation parameter than the BELR, it cannot be part of 
a continuous gas distribution but is rather an independent 
entity. One possibility is that this region is a layer on the 
inner part of the dusty torus that becomes an efficient coro¬ 
nal li ne emitter on ly when evaporated in an X-ray heated 
wind (jPier fe VoitJI 1995l l. We present tentative evidence for 
this scenario in the form of a temperature anti-correlation 
between the coronal line gas and hot dust, which indicates 
that the increased AGN radiation that heats the dusty torus 
appears to increase the cooling efficiency of the coronal line 
gas, most likely due to a stronger adiabatic expansion. 

In a future paper, we plan to present a similar study 
of the coronal line variability for the well-known AGN 
NGC 5548. 
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